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In this study, we demonstrated experimentally and theoretically that the charge transport mechanism in 
amorphous Hf 0 5 Zr 0 5 0 2 is phonon-assisted tunneling between traps like in Hf0 2 and Zr0 2 . The thermal trap 
energy of 1.25 eV and optical trap energy of 2.5 eV in Hf 0 5 Zr 0 5 0 2 were determined based on comparison of 
experimental data on transport with different theories of charge transfer in dielectrics. A hypothesis that 
oxygen vacancies are responsible for the charge transport in Hf 0 5 Zr 0 5 0 2 was discussed. 
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Knowledge about charge transport of high-K di¬ 
electrics is very important for modern microelectron¬ 
ics. Previous transport studies were based on binary 
compound like Hf0 2 and Zr0 2 1-3 . It was shown, 
that charge transport mechanism in Hf0 2 and Zr0 2 is 
phonon-assisted tunneling between traps 2,4 . The charge 
transport mechanism of ternary high-/« solid solution 
Hf 0 5 Zr 0 5 0 2 still remains unknown. 

In this letter, we investigate the charge transport 
mechanism in Hf 0 5 Zr 0 5 0 2 by comparison of experimen¬ 
tal data with different theories of charge transfer in di¬ 
electrics. 

Transport measurements were performed for struc¬ 
tures Si/Hf 0 5 Zr 0 5 0 2 /Ni. To fabricate these structures, 
we deposited the 20-nm-thick Hf 0 5 Zr 0 5 O 2 solid solu¬ 
tion films on n- and p-type Si (10 0) wafers by physi¬ 
cal vapor deposition (PVD). Pure Hf0 2 and Zr0 2 tar¬ 
gets were bombarded by electron beams in high vacuum 
chamber, and Hf x Zr y O were deposited on the wafer form¬ 
ing Hf 0 5 Zr 0 5 0 2 films. A Zr/Hf ratio of ~ 1 was used. 

We did not apply any post-deposition annealing to pro¬ 
duce the most non-stoichiometric films. The structural 
properties of grown high-K Hf 0 5 Zr 0 5 0 2 dielectric were 
examined by grazing incidence x-ray diffraction diffrac- 
togram (GI-XRD). The structural analysis showed that 
the resulting Hf 0 5 Zr 0 5 0 2 films were amorphous. All 
samples for transport measurements were equipped with 
round 50-nm-thick Ni gates with a radius of 70 /im. The 
measurements were performed using a Hewlett Packard 
4155B semiconductor parameter analyzer and an Agilent 
E4980A precision LCR meter. 

Optical (dynamic) permittivity of Hf 0 5 Zr 0 5 0 2 
was calculated in the framework of density functional 
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Figure 1. Experimental (characters) and simulations (solid 
lines) current-voltage characteristics in n-Si/Hf 0 5 Zr 0 5 0 2 /Ni 
structures at different temperatures, (a) Frenkel model (2), 
(b) Hill model of overlapped traps (3), (c) multiphonon trap 
ionization (4), (d) phonon-assisted tunneling between traps 
(5). 


theory using the ab initio simulation code Quantum 
ESPRESSO 5 . Electronic structures and dielectric prop¬ 
erties of monoclinic Hf 0 5 Zr 0 5 0 2 using 12-atom cell was 
simulated. The structure was obtained by replacement 
of a half hafnium atom to zirconium in monoclinic prim¬ 
itive cell of Hf0 2 with following relaxation. The similar 
approach was described earlier 6 . 

Fig. 1(a) shows a set of experimental current-voltage 
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characteristics ( I-V ) of n-Si/Hf 0 5 Zr 0 5 0 2 /Ni structures 
measured at different temperatures T by characters in 
various shapes and colors. Positive applied voltage cor¬ 
responds to positive bias on the Ni contact. The leakage 
current through Hf 0 5 Zr 0 5 0 2 grows exponentially with 
increasing of electric field (or applied voltage) and tem¬ 
perature in accumulation mode (V >0). 

Experiment results were analyzed by using different 
models of charge transport in dielectrics: 

I = eSN 2 / 3 P : (1) 

where / is the full current through the sample, e is the 
elementary charge, S = 7r(70^m) 2 is the contact square, 
N is the bulk trap density, and P is the probability rate 
of charge carrier transfer between traps, which depends 
on the transport model. A mathematical model of well 
known Frenkel law was introduced in 1938 for isolated 
trap ionization 7 : 


P = v exp 
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where v is the frequency factor which was defined as 
v ~ W/h , W is thermal ionization energy of the trap, 
h = 2nH is the Planck constant, j3p is Frenkel coeffi¬ 
cient, F = V/d is the electric field, d is the dielectric 
film thickness, k is the Boltzmann constant, Eoo is dy¬ 
namic permittivity of the dielectric film, and £q is vacuum 
permittivity (electric constant). Results of simulations 
(l)+(2) are shown in Fig. 1(a) by solid lines. One can see, 
that Frenkel model describes the experiment data qual¬ 
itatively very good. However, quantitative fitting pro¬ 
cedure returns underestimated fitting parameter values: 
the slopes of the fitting lines with Frenkel coefficient give 
the dynamic permittivity £oo = 3.6 4- 3.7, which is lower 
than £ oo (Hf0 2 ) = 4.4 8 , £ oo (Zr0 2 ) = 5.6 9 , and calculated 
from the first principals £ oo (Hf 0 5 Zr 0 5 0 2 ) = 4.8 -4 5.2. 
Further fittings return N = 10 7 cm -3 and W = 0.8 eV. 
Found values the charge trap density of N = 10 7 cm -3 at 
v ~ W/h ~ 2 x 10 14 s _1 corresponds to mean distance be¬ 
tween traps s = TV" 1 / 3 ~ 50 /iin that is comparable to Ni 
gate size. Taking all these facts into account one can con¬ 
clude that there is no quantitative agreement between ex¬ 
periments and Frenkel model, despite that Frenkel model 
describes the experiment data qualitatively. 

Simulating in terms of overlapped traps ionization 
(Hill model) 10 


P = v exp 
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are in good quantitative agreement with experiments as 
well as Frenkel model (Fig. 1(b)). However, too low value 
of frequency factor of v ~ 10 7 s _1 was obtained. All val¬ 
ues of obtained filling parameters are collected in Table I. 

Results of simulations by the model of multiphonon 



Figure 2. Configuration coordination energy diagram of trap 
ionization process on negative charged trap in Hf 05 Zr 05 O 2 . 
Lower term is filled ground state, upper term is excited empty 
state. 
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are shown in Fig. 1(c) by solid lines. Here Wph is phonon 
energy, W op t is optical energy of the trap, W t is ther¬ 
mal trap energy, I n are modified Bessel functions, m* is 
the effective mass, and Pi(W) is probability of tunneling 
trough a triangle barrier of W height. Calculated set of 
I-V-T curves is very close to experimental data, obtained 
values of fitting parameters include low trap density of 
N = 2 x 10 13 cm -3 which corresponds to s = 370 nm. 
This mean distance between traps is much greater than 
the film thickness of 20 nm (Table I). Thus, it can be 
concluded that multiphonon trap ionization does not ad¬ 
equately describe charge transfer in Hf 0 5 Zr 0 5 0 2 . 

To get complete vision on the charge transport in 
Hf 05 Zr 05 O 2 experiment data was simulated basing on 
phonon-assisted tunneling between traps 12 : 
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Results of this procedure are shown in Fig. 1(d). The ex¬ 
periment data were described quantitatively and qualita¬ 
tively with the following values of fitting parameters (Ta¬ 
ble I): N = 3x 10 19 cm” 3 , W t = 1.25eV, W opt = 2.5eV, 
TO*/m e = 0.23 (m e is a free electron mass). Fig. 2 shows 
the configuration diagram of a negatively charged elec¬ 
tron trap. A vertical transition with a value of 2.5 eV 
corresponds to the optical trap excitation, transitions of 
1.25 eV correspond to thermal trap energy. 

The same measurements and simulations were per¬ 
formed for p-Si/Hf 0 5 Zr 0 5 0 2 /Ni structures. Results are 
represented in Fig. 3. One can see that all model de¬ 
scribe experimental curves qualitatively. Calculated val¬ 
ues of fitting parameters are summarized in Table I. 
Models of isolated (2) and overlapped (3) charged traps 
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Figure 3. Experimental (characters) current-voltage charac¬ 
teristics and simulations (lines) inp-Si/Hf 0 5 Zr 0 5 0 2 /Ni struc¬ 
tures at different temperatures, (a) Frenkel model (2), (b) Hill 
model of overlapped traps (3), (c) multiphonon trap ioniza¬ 
tion (4), (d) phonon-assisted tunneling between traps (5). 


can describe experiments with inadequate parameters 
like for n-Si-based samples. Multiphonon ionization of 
neutral trap (4) has good agreement with experiments 
at A = 2 x 10 13 cm -3 , which is equal than one got 
for n-Si/Hf 0 5 Zr 0 5 0 2 /Ni structure by the same trans¬ 
port model. At the same time calculated curves in terms 
of phonon-assisted tunneling between traps (5) are close 
to experimental data with the same parameter values as 
that obtained for n-Si/Hf 0 5 Zr 0 5 0 2 /Ni samples. 

Phonon-assisted tunneling between traps adequately 
describes charge transport in Hf 0 5 Zr 0 5 0 2 films on n-Si 
and p -Si substrates. Taking these into account, we con¬ 
clude that the model of phonon-assisted tunneling be¬ 
tween traps describe charge transport in Hf 0 5 Zr 0 5 0 2 
films. Energy parameters of traps, such as the ther¬ 
mal trap energy of 1.25 eV and the optical trap energy of 
2.5 eV, are similar to that in binary oxides Hf0 2 2,13 , and 
Zr0 2 14 . 

Capacitance-voltage ( C-V ) measurements (Fig. 4) 
show that increasing of voltage amplitude leads to shift 
of the hysteresis to negative voltages. This phenomenon 
might be caused by holes trapping on Hf 0 5 Zr 0 5 0 2 in¬ 
terfaces or in the bulk of the dielectric. C-V shift al¬ 
lows us to valuate the density of the filled hole traps as 
n m athrmh h < 3x 10 18 cm~ 3 in the bulk or n m athrmh s < 
6 x 10 12 cm' 2 on the surface states. It is not possible to 
separate percentage of the charge on the surface and in 
the bulk, unusual C-V behavior of reducing the maxi¬ 


Figure 4. Experimental capacitance-voltage characteristics in 
p-Si/Hf 0 5 Zr 0 5 0 2 /Ni structures at different voltage limits. 


mum capacity with an increase in the voltage amplitude 
indicates that the the surface charge is significant. How¬ 
ever, it should be noted that the number of filled bulk 
traps n m athrmh h < 3 x 10 18 cm -3 is much lower than 
total trap density of A = 3 x 10 19 cm' 3 . The possible ex¬ 
planation of this difference is the Coulomb repulsion of 
the charged particles with forming of Wigner-glass-like 
structures 15 in the bulk of Hf 0 5 Zr 0 5 0 2 . 

It was shown that oxygen vacancies are responsible 
for charge transport via Hf0 2 and Zr0 2 2,13,14 . Thus, 
one can expect that oxygen vacancies are responsible 
for charge transport in Hf 0 5 Zr 05 O 2 too. To confirm 
this hypothesis experiments on photoluminescence and 
quantum-chemical simulations are required. 

Recently, it was reported that orthorhombic crys¬ 
talline phase of high-K Hf 0 5 Zr 0 5 0 2 thin films can be fer- 
roelectrics, being perspective material for application to 
ferroelectric random access memory (FeRAM) 16 ' 18 . De¬ 
spite that FeRAM has many advantages, retention char¬ 
acteristics of FeRAM devices much to be desired because 
of the depolarization effect 19 . A possible reason of the 
depolarization effect is charge leakage via traps of the di¬ 
electric. To confirm or refute this hypothesis, transport 
properties of Hf 0 5 Zr 0 5 0 2 in amorphous and ferroelectric 
phases must be compared. 

To summarize, we examined the transport mecha¬ 
nisms of amorphous solid solution Hf 05 Zr 05 O 2 . It was 
demonstrated that all charge transport models such as 
Frenkel model, Hill model, multiphonon trap ionization, 
and phonon-assisted tunneling between traps describe ex¬ 
periment data formally, qualitatively, while only phonon- 
assisted tunneling between traps describes the charge 
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Model 

N (cm -3 ) 

S 

W (eV) 

Wt(eV) Wopt(eV) 

* (a" 1 ) 

£ oo 

m* /m e 

F 

1 x 10 7 

50 /im 

0.8 

— 

— 

~ 10 14 

3.6 4 3.7 

— 

Hill 

3 x 10 19 

3.2 nm 

0.9 4 1.0 

— 

— 

~ 10 6 4 10 7 

5 

— 

MPTI 

2 x 10 13 

370 nm 

— 

0.8 

1.6 

— 

— 

0.17 

PAT 

3 x 10 19 

3.2 nm 

— 

1.25 

2.5 

— 

— 

0.23 

Exp 

~ 10 18 4-10 21 

1 4 10 nm 

~ 1 

~ 1 

~143 - 

- 10 14 4 10 15 

4.8 4 5.2 



Table I. Summary table of the values of the fitting parameters obtained from the simulation I-V characteristics for n- and 
p-Si/Zr0 2 /Ni structures in different models: (F) Frenkel model (2), (Hill) Hill model (Pool law) (3), (MPTI) multiphonon trap 
ionization (4), (PAT) phonon-assisted tunneling between traps (5). The last column represent ranges of expected (reasonable) 
values (from calculation or literature) if any. 


transport in Hfo 5 ZrQ 5 0 2 quantitatively. Comparing ex¬ 
perimental current-voltage characteristics with results of 
simulations revealed energy parameters of the charge 
traps in Hf 0 5 Zr 0 5 0 2 : the thermal trap energy of 1.25 eV 
and the optical trap energy of 2.5 eV, that are equal to 
that for Hf0 2 and Zr0 2 . 

This work was particularly supported by Ministry 
of Science and Technology, Taiwan (grant No. NSC103- 
2923-E-009-002-MY3) (growing test structures, prepar¬ 
ing samples, performing transport measurements), and 
by the Russian Science Foundation (grant No. 14-19- 
00192) (calculations, modeling). 


REFERENCES 

1 Y. N. Novikov, Journal of Applied Physics 113 , 024109 (2013). 

2 D. R. Islamov, V. A. Gritsenko, C. H. Cheng, and A. Chin, Ap¬ 
plied Physics Letters 105 , 222901 (2014), arXiv: 1409.6887 [cond- 
mat.mtrl-sci], 

3 G. Jegert, A. Kersch, W. Weinreich, and P. Lugli, Journal of 
Applied Physics 109 , 014504 (2011), 10.1063/1.3531538. 

4 D. R. Islamov, V. A. Gritsenko, C. H. Cheng, and A. Chin, 
(unpublished). 

5 P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, 
C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococ- 
cioni, I. Dabo, A. Dal Corso, S. de Gironcoli, S. Fabris, 
G. Fratesi, R. Gebauer, U. Gerstmann, C. Gougoussis, A. Kokalj, 
M. Lazzeri, L. Martin-Samos, N. Marzari, F. Mauri, R. Maz- 
zarello, S. Paolini, A. Pasquarello, L. Paulatto, C. Sbraccia, 


S. Scandolo, G. Sclauzero, A. P. Seitsonen, A. Smogunov, 
P. Umari, and R. M. Wentzcovitch, Journal of Physics: Con¬ 
densed Matter 21 , 395502 (2009). 

6 G. Dutta, Journal of Applied Physics 105 , 106103 (2009), 
10.1063/1.3117829. 

7 J. Frenkel, Physical Review 54 , 647 (1938). 

8 S. Kar, ed., High Permittivity Gate Dielectric Materials, 
Springer Series in Advanced Microelectronics, Vol. 43 (Springer, 
2013). 

9 G.-M. Rignanese, F. Detraux, X. Gonze, and A. Pasquarello, 
Physical Review B 64 , 134301 (2001). 

10 R. M. Hill, Philosophical Magazine 23, 59 (1971). 

11 S. Makram-Ebeid and M. Lannoo, Physical Review B 25, 6406 
(1982). 

12 K. A. Nasyrov and V. A. Gritsenko, Journal of Applied Physics 
109 , 093705 (2011). 

13 T. V. Perevalov, V. S. Aliev, V. A. Gritsenko, A. A. Saraev, V. V. 
Kaichev, E. V. Ivanova, and M. V. Zamoryanskaya, Applied 
Physics Letters 104 , 071904 (2014), 10.1063/1.4865259. 

14 T. V. Perevalov, D. V. Gulyaev, V. S. Aliev, K. S. Zhuravlev, 
V. A. Gritsenko, and A. P. Yelisseyev, Journal of Applied Physics 
116 , 244109 (2014). 

15 S. S. Shaimeev, V. A. Gritsenko, and H. Wong, Applied Physics 
Letters 96 , 263510 (2010), 10.1063/1.3458832. 

16 J. Muller, T. S. Boscke, D. Brauhaus, U. Schroder, U. Bottger, 
J. Sundqvist, P. Kucher, T. Mikolajick, and L. Frey, Applied 
Physics Letters 99 , 112901 (2011). 

17 M. Hyuk Park, H. Joon Kim, Y. Jin Kim, W. Lee, 
H. Kyeom Kim, and C. Seong Hwang, Applied Physics Letters 
102 , 112914 (2013). 

18 M. Hyuk Park, H. Joon Kim, Y. Jin Kim, W. Lee, T. Moon, and 
C. Seong Hwang, Applied Physics Letters 102 , 242905 (2013). 

19 C.-H. Cheng and A. Chin, IEEE Electron Device Lett. 35, 138 
(2014). 









